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Abstract 



In the supersymmetric left-right model, the light neutrino masses are given 
by the Type-II seesaw mechanism. A duality property about this mecha- 
nism indicates that there exist eight possible Higgs triplet Yukawa couplings 
which result in the same neutrino mass matrix. In this paper, We work out 
the one-loop renormalization group equations for the effective neutrino mass 
matrix in the supersymmetric left-right model. The stability of the Type-II 
seesaw scenario is briefly discussed. We also study the lepton-flavor-violating 
processes ( r — > /ry and r — > e^ ) by using the reconstructed Higgs triplet 
Yukawa couplings. 



*E-mail: chaowei@mail.ihep.ac.cn 



1 



I. INTRODUCTION 



Current solar, atmospheric, reactor and accelerator neutrino oscillation experiments have 
provided us with very convincing evidence that neutrinos have non-vanishing masses and 
lepton flavors are mixed [1-5]. A global analysis of current experimental data yields 30° < 
012 < 38°, 36° < 9 23 < 54° and < 9 13 < 10° as well as Am 2 21 = m 2 2 - mj = (7.2 • • • 8.9) x 
1(T 5 eV 2 and Am 2 32 = m\-m\ = ±(2.1 • • • 3.1) x 1(T 3 eV 2 at the 99% confidence level 
[6], but three CP- violating phases (i.e., the Dirac phase 5 and the Majorana phases p and 
a ) are entirely unrestricted. These important results indicate that there should be a more 
fundamental theory beyond the standard model, in which three neutrinos are massless Weyl 
particles. One possible candidate for such a theory is the supersymmetric version of the left- 
right symmetric model [7], which provides a natural embedding of the seesaw mechanism 
for small neutrino masses [8]. 

The supersymmetric left-right model [9,10] is based on the gauge group SU(3) C x 
SU(2) L x SU(2) R x U(1) B -l- The quarks and leptons transform under the gauge group as 
Q(3, 2, 1, 1/3), Q c (3*, 1, 2, -1/3), L(l, 2, 1, -1) and L c (l, 1, 2, 1). In the gauge sector, there 
are triplet gauge bosons (W + ,W~ ,W°) L , (W + ,W~,W°) R corresponding to SU(2) L and 
SU(2) R and a vector boson V corresponding to U(1) B _ L , together with their superpart- 
ners. Fermion masses arise from the Yukawa coupling between quarks, leptons and Higgs 
bi-doublets: $ u (2,2,0) and $ d (2,2,0). The gauge group SU(2) R x U(1) B _ L is broken to 
the hypercharge symmetry U(1) Y by the vacuum expectation value (vev) of a B — L = —2 
Higgs triplet A c (l, 1, 3, —2) which is accompanied by a left-handed Higgs triplet A(l, 3, 1, 2). 
The choice of the triplets is preferred because with this choice the seesaw arises from purely 
renormalizable interactions. In addition to A and A c , the model must contain their con- 
jugate fields A and A c to insure the cancellation of the anomalies that would otherwise 
occur in the fermionic sector. Given their strange quantum numbers, the A and A c do not 
couple to any of the particles in the theory, and thus their contributions are negligible for 
any phenomenological studies. The gauge invariant part of the matter superpotential can 
be written as 

W = Y^Q C ) T ^Q + F/(L c ) T i>,L + i [t L T t 2 AL + T c L cT r 2 A c L c ) , (1) 

where $j = ir 2 ^ i is defined and i — u,d. All the couplings Y*, Yf, T and T c are complex 
with T and T c being symmetric matrices. The left-right symmetry implies Y£ = (Y*y 
(a = q, I) and T = T* c . Given the vevs of 0> u ,di A and A c , 

s). w-(s : 



v L oy \v R o. 

the gauge group is broken to U(l) em and the up-type quark, down- type quark, charged lepton 

/ q n ui M d = Yq 



and Dirac neutrino mass matrices turn out to be: M u = Yl L n u , M d = Y / fn d , M l = Y^n d and 



M D = Y l d n u . Meanwhile, the left- and right-handed Majorana neutrino mass matrices can 
be obtained from the corresponding mass terms in Eq. (1) once the Higgs triplets A and A c 
acquire their vevs: M L ~ v L T and M R ~ v R T . Integrating out the heavy particles (i.e., the 
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right-handed Majorana neutrinos and Higgs triplet), one obtains the effective mass matrix 
for three light (left-handed) Majorana neutrinos via the Type-II seesaw mechanism [11]: 

M v ~ M L - M T D M^M D ~ v L T - —MlT~ x M D . (3) 

V R 

We may find that the same coupling T appears in both contributions just because of the 
left-right symmetry. 

Note that Eq. (3) has a duality property [12]: given M D , there exist eight possible Higgs 
triplet Yukawa couplings which result in the same neutrino mass matrix. The stability of the 
duality relation and some other phenomena based on this have been investigated recently. 
In this paper we perform a full analysis of the renormalization group equations (RGEs) 
of the effective neutrino mass operators. We write down the /3-functions of the effective 
neutrino mass operators and discuss the stability of the Type-II seesaw mechanism. Lepton- 
flavor-violating decays in the supersymmetric left-right model are different from that in 
the minimal supersymmetric standard model (MSSM) for the existence of the Higgs triplet 
Yukawa coupling T [13,25]. In this paper, we calculate the BR(r — > fxrf and BR(r — > erf 
by using the reconstructed Higgs triplet Yukawa couplings in the supersymmetric left-right 
model. 

The remaining part of this paper is organized as follows. In section II, we calculate the 
one-loop RGEs for the effective neutrino mass operators. Section III is devoted to studying 
the lepton-flavor-violating processes. A summary of our main results is given in section IV. 
Some useful formulas are listed in appendices A and B. 



II. RENORMALIZATION GROUP EQUATIONS OF THE EFFECTIVE 
NEUTRINO MASS OPERATORS 



We assume that the gauge and discrete left-right symmetries are both broken by the vev 
of A c at the high energy scale in our model. As a result the right-handed neutrinos and Higgs 
triplets are much heavier than other particles. Integrating out the right-handed neutrinos in 
the leading-order approximation, one obtains the effective neutrino mass operators, which 
are contained in the F-term of the superpotential, 
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-M,)Jle c ^ u )J f a e a \<S> d ) bl + h. C ., 



(4) 



where 



k 2 — 2 



(5) 



Due to the non- renormalization theorem [14], the RGEs for operators of the superpotential 
are governed by the wave function renormalization for the superfields. At the one-loop level 
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the wave-function renormalizaton constants Z are obtained with the dimensional regular- 
ization via the dimensional reduction [15]: 



6Tr 
6Tr 
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{Y-)X] + 2Tr [(^)t V] - ^ 2 - 3^ 2 2 , 
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(6) 



Using the counterterms calculated above and the technique described in Ref. [16], we obtain 
the /9-functions {(3 X = H>-$-X) of the effective mass operators n i [i = 1,2,3) and the Higgs 
triplet Yukawa coupling T : 



16tt 2 /3 Ki 



16tt 2 # 



n T ■ «! + «! • ft + j 6Tr [(V) f y; 

ft T • K, + /,', • ft - <J I, ! " 

ft r - 
ft T 



WW 

K 3 + K 3 .n+{3Tr (Y q u )% u 
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5^ 
6 2 

3Tr 
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+ J 7 • ft + 2Tr 



9 2 
5* 



6 2 

5^ 



6# 
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where 



ft = ( I7)ty« + {Y d )% d + 

Some comments are in order. 



(7) 



(8) 



• In calculating the /3-functions, we have assumed M A (the mass of the Higgs triplet) to 
be lighter than Mi which is the mass of the lightest right-handed neutrinos. Actually 
this assumption is not necessary. One may integrate out u R and A each at its own mass 
scale and redefining iteratively the effective operator, which is more reasonable. Below 
m A , the /9-functions of the effective mass operators, which come from integrating out 
the Higgs triplet, are similar to ft/s. 

• Given the vacuum expectation values of the Higgs bi-doublets and triplets in Eq. 
(2), only Ki gives rise to masses of the light left-handed neutrinos after spontaneous 
electro- weak symmetry breaking. We just need to calculate the /3-function of k 1 when 
considering the renormalization group effects of neutrino mass operators. Besides, all 
operators in W K contribute to the lepton-flavor-violating processes. However, such 
processes are strongly suppressed by heavy masses of the right-handed neutrinos. 

• Below the lightest seesaw scale, the /9-function of the effective neutrino mass operator 
prossess the same as that of the Type-I seesaw model in the MSSM, only up to a 
replacement Y?Y t — ► {Y t u )% u + {Y d )% d . 
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Due to the renormalization group (RG) evolution effects between the M A and M l scales, 
the seesaw formula in Eq. (3) is modified, where two JF's in Type-I and Type-II terms are 
not equal anymore. As a result the duality property is slightly broken when considering the 
RG evolution effects of T and the effective neutrino mass operator. 



III. LEPTON FLAVOR VIOLATION IN THE SUPERSYMMETRIC LEFT-RIGHT 

MODEL 

In this section, we first give the analytical formulas to be used for the calculation of the 
lepton-flavor-violating processes and then list our numerical results. 



A. Analytical formulas 



Working in the basis where the sleptons are in weak eigenstates together with the 
charginos (neutralinos) in their mass eigenstates, we write down the interaction Lagrangian 
of lepton-slepton-chargino in the following form: 



^int — ^ V \iX.A^LR 



+h.c. 



iMi) td I r^W p \; I .~.t in A( S) p _i_ n A ® 

LR r R ' ^LL r L) L i ' v RiXAV^RR r R + u RL 

^PR + Nti ] PL)h + ~ei^ 



RiXAV^RR r P + JV Rf. 



RL ^L)h 

t rTPl% 



(9) 



where the coefficients are: 
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(10) 



Here O l , O r and O n are real orthogonal matrices that diagonalize chargino and neutralino 
mass matrices respectively. Their explicit forms are listed in appendix A. tan/3 = kJ K d is 
defined. 

Let us discuss the branching ratios of the lepton-flavor-violating processes in the super- 
symmetric left-right model. The radiative decays l { — > Z • + 7 are induced by the effective 
operator [17]: 
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el 3 {iD\P L + iD\P R ) an.F^ + h.c. , (11) 

where e and F are the charge and the electromagnetic field strength, respectively. These 
operators are chirality-flipping (dipole) and come from SU(2) h x U(\) Y - invariant operators 
with at least one Higgs field. 

In the "mass insertion" method and leading-log approximations, the coefficients D\ R can 
be calculated [13] and we write down the explicit expression in appendix A. The branching 
ratio of U — > lj + 7 decay due to the new contributions is given by: 

BR(/ l - Z i7 ) = ^ (I Dl | 2 + I Dl | 2 ) BR& — Z^-i/J, (12) 

where a = e 2 /(47r), G F is the Fermi constant, BR(r — > Atz/ T z7 M ) 17% and BR(r — > ez/ T 77 e ) « 
18% [18]. 

In the minimal SUGRA scenario, at the gravitational scale the supersymmetry breaking 
masses for sleptons, squarks and the Higgs bosons are universal, and the SUSY breaking 
parameters associated with the supersymmetric Yukawa couplings or masses are proportional 
to the Yukawa coupling constants or masses. Then, the SUSY breaking parameters are given 
as: 

( m l)ij = ( m R)ij = ( m l)ij = $ij m l, 

m l, = m l 2 = m l, 

(A^ d r = (Y l u ' d ra ,A i ^ = ^a , 

B'J = M ViV .bo,B*=fibo. (13) 

Flavor violation in the slepton sector arises from radiative corrections induced by the flavor- 
violating couplings of heavy states populating the theory between the Planck scale and the 
electroweak scale. Integrating the one-loop renormalization group equations [19] for the soft 
breaking masses m|, m R and trilinear A^' d in the lowest-order approximation, one obtains 
the off-diagonal term for m|, m R and A^ d : 

(mlh - (m|),, * -^fj^' - -i^'V^' (14) 

where 

K, = [lTW)* + Y!W)% log + 3(^%log . 

These off-diagonal terms generate new contributions in the amplitudes of lepton-flavor- 
violating processes [20] such as r — > /I'j and r — > e7 . 



B. Numerical results 

The lepton flavor mixing matrix (?7 MNS ) comes from the mismatch between the diagonal- 
izations of the neutrino mass matrix and the charged lepton mass matrix. The tri-bimaximal 
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mixing pattern [21] is strongly favored by the solar and atmospheric neutrino oscillation 
measurements: 



/ -2- 



MNS 



1 

1 1 



m \ 
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V y/E 



1 

^3 



1 

1 

y/2/ 



(15) 



A global analysis of current experimental data yields the values for the solar mass splitting 
Am 2 12 = (8.0±0.3)xlO~ 5 eV 2 and the atmospheric mass splitting \m 2 23 \ = (2.5±0.2) x 10~ 3 eV 2 
[6]. We assume that three light left-handed Majorana neutrinos are in normal mass hierarchy 

so that rrin z 



Am| 3 



0.05eV and m 2 ~ 



Am 



12 



0.009eV. We 



(i.e., m 1 < m 2 < m 3 

also take m 1 ~ O.OOleV, v L ~ 0.05eV and v L v R /vf l ~ 1, which are natural values 

We assume that at the GUT scale the theory is given by the supersymmetric SO (10) 
model which contains two 10-dimensional and a pair of 1260126 representation Higgs bosons. 
Then the most general Yukawa couplings lead to the following mass relation for the fermions: 
M u = M D . We neglect the CKM relations between the up- and down-type quarks in our 
numerical calculations, assuming that the up-type and down-type quark mass matrices are 
both diagonal. The Dirac neutrino mass matrix turns out to be M D = diag(m u , m c , m t ). 

Using these choices and the technique described in [12], one obtains eight different solu- 
tions for the triplet Yukawa coupling T through the left-right seesaw formula in Eq. (3): 



/ -0.00169 -0.00349 0.00015 \ 
-0.00349 0.51022 -0.51309 

V 0.00015 -0.51309 0.69097 / 

/ 0.06235 0.06316 0.05950 \ 
0.06316 0.04996 0.07515 , 

\0.05950 0.07515 0.21616/ 

-4-10- 10 
4 • 10" 8 
6 • 10- 6 

5- 10- 11 
-3- 10~ 8 
-6- 10~ 6 




( 0.06236 
0.06316 
V 0.05952 

-0.00169 
-0.00349 
0.00016 

0.06067 
0.05967 
0.05966 

/ 0.06067 
0.05966 
\ 0.05967 



0.06316 
0.04995 
0.07326 

-0.00349 
0.51021 
-0.51498 

0.05967 
0.56017 
-0.43888 

0.05966 
0.56016 
-0.44078 



0.05952 \ 
0.07326 
-0.13080 / 

0.00016 
-0.51498 

0.34400 

0.05967 
-0.43888 
0.73374 

0.05967 
-0.44078 
0.38678 



(16) 



It is easy to check that the duality relation (J 7 i +J 7 ' i — m u /v L ) is satisfied very accurately for 
the solutions given above. 

Now, we present our numerical results of BR(r — > //, e + 7) in the parameter space given 
above. The experimental upper limits on those branching ratios are: BR(r — > /i + 7) < 6.8 x 
10~ 8 and BR(r -> e + 7) < 1.1 x 10~ 7 at 90% C.L. [23] and the sensitivities of a few planned 
experiments [24] may reach BR(r -> e + 7) ~ cjlO" 8 ) and BR(r -> fj, + 7) ~ C(10- 8 ). 
FIG. 1 and FIG. 2 show the BR(r — > [/x, e] + 7) changing with m . We find that the 
experimentally allowed ranges of BR(r — > [//, e] + 7) can be reproduced from all of these 
eight different triplet Yukawa couplings in the chosen parameter space. Besides, curves 
corresponding to T z and F A are lapped over with each other because there is little difference 
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in their numerical expression. Although eight different Higgs triplet Yukawa couplings result 
in the same neutrino mass matrix through the Type-II seesaw formula, their effects on lepton- 
flavor- violating processes are very different. As a result, we may check the stability of the 
Type-II seesaw formula by measuring the branching ratios of the lepton-flavor-violating r 
decays accurately in the future experiments. 

IV. SUMMARY 

In addition to the right-handed neutrinos, the Higgs triplet is another source of the 
neutrino mass generation in the Type-II seesaw model, so the evolution of the neutrino 
mass matrix is a little different from that in the Type-I seesaw model. Besides, the duality 
property for the Type-II seesaw formula indicates that there exist eight possible Higgs triplet 
Yukawa couplings T which , for a given M D , result in exactly the same mass matrix of light 
neutrinos. In this article, we have calculated the RGEs for the evolutions of the Type- 
II seesaw neutrino mass matrices from the seesaw scale to the electro-weak scale in the 
supersymmetric left-right model. Instead of giving numerical analysis, we have discussed the 
stability of the Type-II seesaw model. On the other hand, the Higgs triplet Yukawa coupling 
is an important source for the lepton-flavor-violating r decays. We have calculated these 
eight Yukawa couplings through the Type-II seesaw formula and applied them to evaluating 
the branching ratios of lepton-flavor-violating r decays. We find that their contributions to 
the branching ratios are different and the stability of the Type-II seesaw can be checked by 
measuring rare r decay accurately. 

In conclusion, the supersymmetric left-right model supplies an interesting platform for 
the neutrino sector, which could be tested in the future LHC and ILC experiments. 

ACKNOWLEDGMENTS 

The author is indebted to Professor Zhi-zhong Xing for reading the manuscript with 
great care and patience, and also for his valuable comments and numerous corrections. He 
is also grateful to S. Zhou and H. Zhang for useful discussions. This work was supported in 
part by the National Nature Science Foundation of China. 



8 



APPENDIX A 



In this appendix, we consider chargino mixing and neutralino mixing in the supersym- 
metric left-right model. We first write down the A — — A terms of the Lagrangian, which 
involve the soft supersymmetry-breaking terms and the scalar potential [9,26]. 

£ GH = +iV2Tr[(a ■ A L )^{g L a ■ X L + 2g v X v )a ■ A L ] + h.c. 
+i v / 2Tr[(a • A R )\g R a ■ X R + 2g v X v )a ■ A R ] + h.c. 

+ -^=Tr[$t (g L a ■ X L + g R a ■ A R )$ M ] + h.c. 

+^=Tr[^ d (g L a ■ X L + g R a ■ X R )$ d } + h.c. 

+Tr[fi 2 (a ■ A L )(a ■ 5 L )} + Tr[^3(^ • A R )(a ■ 5 R )} + h.c. 
+m L (X a L X a L + X a L X a L ) + m R (X a R X R + X R X R ) + m v (X v X v + X V X V ) 
+Tr[ f i 1 (a 1 $ u a 1 ) T $ d }. (Al) 

Substituting the vacuum expectation values of the Higgs fields from Eq. (2) into Eq. (17), 
Keeping only the terms involving charged fields, we get 



?c = {^ R (V2g R v R A R + g R k d 4>\) + iX L (V2g L v L A[ + g L k d 

+iX R g R k u (f)- + iX{g L k u (f)- + 4m L X{Xl + Am R X R X R + fn^ d 
+Vi<f>Z4>d + ^At~5 L + fi 3 A R 5 R ] + h.c. 



(A2) 



We consider the chargino mass matrix M c , which is a 6 x 6 matrix appearing in the chargino 
mass terms. 



In this model, ip is defined to stand for the following fields: 

V + = (-^,-a+,0+ 0+A+A+) r , 

ip- = (-i\£, -iX R , <f)~,(t) d ,5l, 5 R ^j . 
Comparing Eq. (19) with Eq. (17), we write down the explicit expression of M c : 
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A*3 J 








, we 


can diagonalize Mq by 6 x 6 



(A6) 



O r and O l according to O r M c O\ = M§, where Mq is a diagonal matrix. It is tedious 
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to write down the analytical expressions of O l and O r . 
expressions: 



Hence we only list their numerical 
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-0.075 
-0.050 
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-0.001 
-0.034 
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0.005 
-0.008 
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0.009 


-0.739 
-0.644 




-0.002 
0.090 


0.817 
0.570 







-0.371 



0.929 
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-1 




0.001 
0.062 


0.656 
-0.752 




-0.009 \ 
-0.002 


-0.572 
0.820 



0.001 
0.929 





0.371 



0.981 \ 
-0.002 


0.147 
0.129 
/ 



(A7) 



Here we choose M L = 1 TeV, M R = 20 TeV, [ii = /i2 



H3 



200 GeV, tan/? = 1.5, 



0.05 eV and v 



R 



10 GeV in our calculation. 



In order to obtain the neutralino part of the Lagrangian, we replace the vevs of the Higgs 
bosons into Eq. (17), keeping only the neutral terms: 



V -N 



{-iV2(X° L g L - 2X° v g v )v L A° L - i\ f 2(X R g R - 2X° v g v )v R A R 

+i ^ X ° R9R ~~ x °l9l)^u(Pi u - i-^{\ R g R - X° L g L )K d ^° 2d 
+m L (X° L X° L + A° L A°) + m R (X R X R + A^) + m v (X° v X° v + X° V X° V ) 



-/•til 



&&)} 



h.c. 



(A* 



The neutralino particles are produced in two stages of symmetry breaking [27]. The first 
stage, the vev v R is responsible for giving masses to the heavy neutralinos. The second stage, 
the vevs k u and K d are responsible for giving masses to the light neutralinos. The amount 
of mixing between heavy and light neutralinos is small, so one can calculate the neutralino 
mass eigenstates for both stages as independent cases. 

We define £ N : 

U = HA£, -ix% 4>° lu , 4 U , 4 d , 4>° 2d ), (A9) 

Then the relevant part in Eq. (24) may be written as: 



(A10) 
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M N is diagonalized by a real orthogonal matrix O n with (D N M N Ojj = M$. We write down 
the numerical expression for O n : 



( 



o N = 





-0.995 
-0.106 





0.018 



Here we choose M L = 
calculation. 



0.999 
0.001 
-0.006 





0.001 



0.005 
-0.088 
0.707 





-0.702 







-0.707 
-0.707 








0.707 
-0.707 




0.004 \ 
-0.061 
0.700 





0.711 ) 



(A12) 



1 TeV, M R = 20 TeV, ^ = 200 GeV and tan/3 = 1.5 in our 



APPENDIX B 



In this appendix, we write down the formula of D\ R 1 , which are a little different from 
the formula given in Ref. [13]: 

n7 1 _ M . oN ^) N A( j ) A i l( 2 V . ( 1 1 g n {Ml /mlJ 

~ 2W M *° NrR NlL K UL - mi u m* ~ m\ Li m\ 



a ( M 2 .. /rh.. 

e L j> 



1 1 ^(Mjo/mlJ 1 1 g n (M? /m 2 



2222 2 '2 22 2 2 

m~ T . — m~ B . m~ T . — m~ T . m~ r . m~ T . — mj D . mf r . — mj r . m~ r . 



>2 e i LL LL -2 -2 -2 " ' ' ~' -2 

Y mf. —m~.\ m~. m~. 

1 -M-oN A(t) N AU) ^ (9n(M 2 ^/ml) g n {M^lm\ 



1 r^C^A'Hm 2 ) ( 1 1 g c (Ml/mlJ 

v ' \ ^fli ^Li v Ri v Lj v R,i 

i l gcfff-My _^ i g.iM; /< 

noon o 



m?. — m?. ml — to?, to?, ml — m?. m?. — m?. m?. 

, 1 m ^ cm (™lh ( fc(M}-/ml) f c {Ml/ml 3 ) \ 

~Ci7~A \2 e i LL ^LL ~~2 ~ T ~ 2 ~ 2 

6(4vr)^ m$.—m$. \ m~. m~. / 

(4tt) 2 x lk "i?-m?. V m?. m\ 



n i _ 1 A/f n N A(i) M A(j) A ii,-2\ ( 1 1 g n (M^/ml R J 



we do not consider the contributions of the double charged chargino mediated diagrams, since 
their contributions are very small. 
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g n (M? /fhlJ + 1 1 g n (M} /ml R .y 



—m N v ' N yj> - eR/ - x 

2 e * RR RR -L -I, y 

1 McN^N*® f 9 "^ 1 ™^ 9n(Ml /mlJ - 



2(4-) 2 * M M fnl m -fnl Rj \ m\ m m 2 Rj J 

I 1 M r A(i) r AU) A u (rh 2 )-- ( 1 1 g c (M?-/mlJ 

1 1 9c(Ml/mlJ 1 1 ft^-MS^ 

noon o 



— m 2 , — ml ml ml — ml ml — ml ml 

v Li V R% v Li v Rj v Li v R,j v Ri v Rj v Li v Rj 

+ 6(4-) 2 L ml m - ml R] \ m\ m m 2 Rj ) 

where the loop functions are 

1 

2(1 ~ x y 



fn(x) = -r(2 + 3x — 6x 2 + + 6rrloga;), 



/ c (a;) = - — 1 . 4 (2 + 3x - 6x 2 + x 3 + 6x\ogx), 
2(1 — xf 

9c{x) = ]_ (3 - Ax + x 2 + 21ogx). (B3) 
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FIGURES 
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FIG. 1. Illustrative plot for BR(r — ► /i+7) changing with m . We take tan (3 = 1.5, M L = 1 Tev 
and M R = 20 Tev in our plot. Here the dot line corresponds to T x \ the dash dot line corresponds 
to J-\; the short dash line corresponds to T 2 \ the short dash dot line corresponds to T 1 2 \ the solid 
line corresponds to J- 3 and the dash line corresponds to the short dot line corresponds to 
T 1 ^ the dot horizontal line corresponds to the experimental upper bound. 




FIG. 2. Illustrative plot for BR(t — ► e+7) changing with m . We take tan (5 = 1.5, M L = 1 Tev 
and M R = 20 Tev in our plot. Here the dash line corresponds to T x \ the dot line corresponds to 
T\\ the dash dot line correspond to J- 2 ] the dash dot dot line corresponds to J-' 2 ', the solid line 
corresponds to the short dash line corresponds to J-\] the short dash line corresponds to 
and J-^; the dot horizontal line corresponds to the experimental upper bound. 
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